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Quark number susceptibilities approach their ideal gas limit at sufficiently high temperatures. As in the 
case of other thermodynamic quantities, this limit itself is altered substantially on lattices with small temporal 
extent, Nt = 4-8, making it thus difficult to check the validity of perturbation theory. Unlike other observables, 
improving susceptibilities or number densities is subject to constraints of current conservation and absence of 
chemical potential (/x) dependent divergences. We construct such an improved number density and susceptibility 
for staggered fermions and show that they approximate the continuum ideal gas limit better on small temporal 
lattices. 
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1. INTRODUCTION 

Thermodynamic observables, such as the en- 
ergy density e, the pressure P or the quark num- 
ber susceptibility x, approach the corresponding 
ideal gas values as the temperature T (or the 
chemical potential /i) becomes very large. Cor- 
rections to the ideal gas behaviour can be com- 
puted in a weak coupling approximation. Such 
a picture of ideal or weakly interacting plasma is 
often used in most physical applications in heavy 
ion phenomenology and in cosmology of the Early 
Universe. Whether the quark-gluon plasma can 
be described by such a simple picture in the rel- 
evant range of temperatures, say 1 < T/T c < 10, 
is a very important issue for these applications, e 
and P are known to show strong deviations from 
the ideal gas picture in this range which cannot 
be explained by naive perturbation theory up to 
0(g 5 ). Clever resummations of the perturbation 
series are thought to be a cure. One such ap- 
proach (T) successfully describes the lattice data 
for P above ~ 3T C . The quark number suscepti- 
bilities provide an independent test of any such 
resummation scheme. Recent lattice results j| 
for these exhibit a strong disagreement with both 
the naive perturbation theory as well as the re- 
summed H one. These lattice results were ob- 
tained using staggered fermions on lattices with 



4, 6 and 8 temporal sites in a temperature range 
1 < T/T c < 3. The ratio x/Xideai from these cal- 
culations were extrapolated to the continuum to 
compare with the perturbative results. 
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Figure 1. The susceptibility x/T 2 as a function 
of N t for 7V S 3 x N t lattices with N s = 3N t . 

As shown in Figure [l], the ideal gas result, 
Xideai, for the naive staggered fermions varies 
strongly for Nt =4, 6 and 8, and deviates sub- 
stantially from the continuum result, indicating 
the need for an improved operator for it in order 
to confirm the implied non-perturbative nature of 
the quark-gluon plasma by the results of || . Re- 
call that x ( or the number density) is obtained 
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from In Z by taking second (or first) derivative 
with respect to fj,. Improving it therefore needs 
to be done by modifying the number operator TV, 
which in turn is fixed by the (baryon) current 
conservation equation. Following the study || 
of improvement of the energy and pressure, we 
employ the Naik action Q and the P4 action Q 
and write down N for them from the correspond- 
ing current conservation equations. An additional 
complication for nonzero fi is the presence of di- 
vergences even for the ideal gas. Prescriptions Q 
were proposed for the naive action to eliminate 
them, which were shown |tJ to satisfy a general 
condition. We demonstrate that the same condi- 
tion can suffice for a class of improved actions as 
well. Numerical results for the ideal gas are pre- 
sented to contrast %(0, T), n(0,T) and n(/z, T) 
with the corresponding continuum results. 

2. IMPROVED NUMBER DENSITY 

We follow the notation of Ref. for the gen- 
eral form of the fermion actions, which yields the 
naive, the Naik and the P4 action for their Ci o, 
c 3 .q and ci )2 equal to (^,0,0), (jq,j$, 0) and (§, 0, 

) respectively, where these coefficients multiply 
the first, third and mixed third derivatives respec- 
tively. Setting all the link variables to unity, the 
current conservation equation for the free theory 
can be easily shown to be : 



c l,2 Ei///j j A M +2i/j^ + ^\i-Zv3$ 

where the currents are defined by 



0, (1) 



1pxJ^x+p,-2i> + 4>x+jl~2vl^x- (2) 



The A's in eq. ([j]) denote appropriate backward 
difference operators. Summing eq. ([j]) over all 
spatial lattice sites x, one obtains ^2^N(x) = 
constant, where the "number density" 

N(x) = c 1>0 j£\x)+3c 3fi jP(x)+ 



2ci, 2 E^obS ) W-iS ) W]+ 

ci,2 E^obo? (*) + iS (*)] + A tF(x) . (3) 

Adding the term /x J2 X N(x) to the /i = action, 
the last term drops out due to the anti-periodic 
boundary conditions on the if} and tp in the time 
direction. Following Ref. , the improved action 
for /i ^ can thus be obtained by introducing 
functions f n (/ia) = 1 + n/ia (g n (fj,a) = 1 — nfia) 
for n =1, 2, and 3, which multiply the terms con- 
taining forward (backward) links in the time di- 
rection. The coefficients in eq. (||) dictate that 
/„ (or g n ) multiplies an n time-link term. 

Without giving the details here, we assert that 
a) the improved action with these f n and g n does 
lead in general to /J? /a 2 {n/a 2 ) divergences in the 
energy (number) density and b) the ansatz /„ = 
/"> ff" = 9i with fi(na) ■ g\ (fia) = 1, eliminates 
them. This ansatz is a natural generalization of 
the results in for the naive action. Thus, as in 
the naive fermion case ||, /i(M a ) — exp(/ia) or 
fi (fia) = (1 + /id) / y/l — fi 2 a 2 can be used for the 
improved action in conjunction with the above 
ansatz for f n and g n . 

3. NUMERICAL RESULTS 

Figure [j] compares the behaviour of the ideal 
gas susceptibility at fj, = 0, Xideah f° r the Naik 
and the P4 action with that for the naive ac- 
tion. The results have been obtained on lattices 
with an aspect ratio of 3 (N s = 3N t ) and for a 
fixed quark mass m/T. Since f'(0) — —g'(0) = 
/"(0) = g"(0) = 1 for all the allowed / and g, the 
results in Figure [l] are valid for all prescriptions 
of /. One sees that both the improved actions 
result in a much milder A t -dependence in Xideai 
and the results are close to continuum by N t — 8 
already. Interestingly, the improved action results 
approach the continuum result from the opposite 
direction as compared to the naive action. 

As one switches on /i, the susceptibilities for the 
improved actions continue to be close to unity on 
lattices with Nt > 8 but they do exhibit a mild 
/i-dependence, increasing above unity for fi *~ T. 
For small /i one also sees a similar pattern in the 
number density, n/T 3 , as seen in Figure |: the 
improved action results are closer to the contin- 
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uum results than the naive action for Nt = 8. 
The results also show a mild /^-dependence which 
may be a cause of concern for the methods which 
attempt to study the equation of state at \i =0 
by using a Taylor expansion in /i in the simula- 
tions and could result in systematic effects in the 
predictions for the plasma phase. 
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Figure 2. Comparison of the continuum number 
density as a function of /i/T with the naive, Naik 
and P4 improved actions on 12 3 x 4 and 24 3 x 8 
lattices. 
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Figure 3. Same as Fig. | but on 24 3 x 48 lattices. 

Figure || shows the results obtained for the 
number density of an ideal gas at very small tem- 
peratures by taking N t very large. On a 24 3 x 48 
lattice, the number density seems to be essentially 
the same for small /i for all actions, including the 
naive action. However, none of them has the cor- 
rect continuum /i-dependence. Only when /x/T 
is large do they all approach the continuum. Im- 
proved actions seem to be closer to the continuum 
result for large /i. Note that even at the highest 
/i/T, /ia is still only about 0.6 in Fig. |[ The 



deviations for small /i/T are due to the smallness 
of the spatial volume. 

4. SUMMARY 

We have constructed an appropriate number 
density term, N, for a set of improved actions 
from the current conservation equations for these 
actions. The Naik action and the P4 action, 
which correspond to specific choices of values for 
the action parameters, are included in this set. 
Adding fj,N to the action for /i = 0, where /i is the 
baryon chemical potential, leads to divergences 
in the energy density and the number density, 
as for the unimproved naive action. These can 
be eliminated by the same condition as for the 
naive action, namely, /(/i) ■ g(fi) = 1, provided 
for each forward (backward) time link the corre- 
sponding gauge variable Uq(x) (£/q(x)) is multi- 
plied by /(/i) (g(p)). 

The quark number susceptibility and the num- 
ber density for the ideal gas do approach the 
continuum result on smaller temporal lattices for 
the improved action. It would be interesting to 
study them in the high temperature regime of 
QCD. Finite size effects appear to be mildly /i- 
dependent for the susceptibility and the number 
density which may lead to systematic effects. 
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